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INTRODUCTION 


Physical  vapor  deposition  produces  structures  that  are  topologically  and  morphologically 
metastable  (ref  1).  Metal  vapors  condense  into  a  fine-grained  crystalline  form  on  a  chilled 
substrate,  resulting  in  evolving  structures  that  are  inherently  more  disordered  and  less  relaxed 
than  those  produced  by  other  coating  processes.  The  process  control  parameters  determine  the 
coarseness,  shape,  and  distribution  of  these  microstructures.  These  microstructural  features 
affect  the  material  properties  and  therefore  the  performance  of  a  coating  (ref  2).  These  features 
can  be  quantified  in  terms  of  self-affine  scaling  parameters  (refs  3,4)  that  provide  details  on  the 
intrinsic  structure  of  the  coating.  Furthermore,  the  time  evolution  of  these  scaling  parameters 
defines  a  unique  universality  class  that  is  associated  with  the  deposition  process,  and  provides 
insight  into  the  dynamics  of  the  growth  processes  of  thick  metal  films. 

APPROACH 

Sputtered  niobium  and  tantalum  coatings  were  prepared  for  analysis  at  various  deposition 
stages  using  different  substrates.  Niobium  was  sputter  deposited  on  silicon  wafers,  while 
tantalum  was  sputter  deposited  on  silicon  and  stainless  steel.  The  evolving  coating  morphology 
was  characterized  in  terms  of  dynamic  scaling  exponents  using  atomic  force  microscopy  (AFM) 
to  map  surface  structures  over  a  range  of  scales  from  10-nm  to  5-fXm.  The  coating  surface  was 
assumed  to  be  a  homogenous,  self-affine  structure  that  is  statistically  invariant  under  anisotropic 
dilations  for  the  AFM  scan  size  employed  in  the  measurements.  Anisotropic  scaling  was 
measured  in  terms  of  the  roughness  exponent  a,  parallel  correlation  length  ^//(t),  and  the 
perpendicular  correlation  length  %±(t).  These  scaling  parameters  provide  information  relating  to 
the  intrinsic  properties  of  the  surfaces  over  a  range  of  scales  after  t  minutes  of  sputtering.  The 
anisotropic  scaling  described  by  ar  applies  over  a  scaling  range  that  is  measured  in  terms  of 
^//(t).  The  parallel  correlation  length  is  the  distance  beyond  which  there  is  no  correlation  in 
heights  between  points  on  the  surface.  The  perpendicular  correlation  length,  %±{t),  characterizes 
fluctuations  in  the  growth  direction.  It  is  related  to  the  RMS  surface  roughness,  cr(t),  by  £jt)  =  V 
2  *  a(t).  The  values  of  a,  £//(t),  and  gjt)  were  determined  using  a  generalized  form  of  the  height 
correlation  function  (ref  4) 


Ch(r, t )  =  ([h(r0+r,  t+t0)  -  h(r0t0)]2  >  ,o,ro  .  ( 1 ) 


where  h(r,t )  is  a  single-valued  height  of  the  surface  at  location  r  at  time  t. 


Correlations  between  points  are  time-dependent  and  generally  increase  with  sputter  time. 
Increasing  grain  size  may  be  characterized  by  growth  in  the  lateral  (£//t))  direction  by  the 
dynamic  exponent  z.  Corresponding  changes  in  the  growth  direction  are  described  by  the  growth 
exponent,  (3.  In  general,  due  to  anisotropy  of  the  deposition  process,  [5  According  to  the 
dynamic  scaling  theory,  at  small  sputter  times,  £//(t),  and  %±(t)  are  given  by 


& 

(2) 

(3) 
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The  scaling  exponents  a ;  /?,  and  z  define  a  unique  universality  class  that  is  independent  of 
the  details  of  the  deposition  process.  They  provide  a  unique  metric  for  describing  the  evolving 
surface  structure  of  the  coating. 

Ch(r,t )  is  sensitive  to  dominant  features  of  the  substrate  and  variations  in  the  curvature  of 
the  sample  (ref  3).  The  influence  of  these  effects  can  be  minimized  by  analyzing  the  structure 
relative  to  a  curved  reference  surface  (refs  5,6)  or  by  dividing  the  surface  into  patches  that  are 
analyzed  separately  using  the  "min-max"  method  (ref  7).  Meisel  (ref  5)  suggests  a  technique  for 
measuring  scaling  properties  relative  to  a  curved  reference  surface  by  limiting  analysis  to  local, 
statistically  homogenous  regions  that  are  large  enough  to  define  a  statistically  reliable  Ch(r,t). 
However,  determining  the  size  of  local  reference  surfaces  that  properly  adapt  to  the 
heterogeneity  of  the  surface  requires  subjective  investigator  input.  The  "min-max"  method  (ref 
7)  requires  the  investigator  to  judge  the  number  of  local  maxima  and  minima  of  the  surface  to 
use  in  the  analysis. 

Therefore,  the  time  evolution  of  the  coating  surface  was  also  characterized  in  terms  of 
scaling  exponents  using  detrended  fluctuation  analysis  (DFA)  (ref  8)  to  correct  for 
inhomogeneities  in  the  structure.  Peng  introduced  DFA  to  determine  long-term  correlations 
embedded  in  time-series  data.  The  sequence  of  integrated  data  of  length  N  is  divided  into  Nil 
nonoverlapping  segments  of  length  l.  The  average  local  variance  of  the  segments  is  measured 
relative  to  a  linear  least-squares  fit  to  the  data  within  each  segment  /.  The  dependence  of  the 
variance  on  /  defines  the  scaling  properties  of  the  time-series.  In  the  present  study,  DFA  was 
used  to  determine  the  scaling  properties  of  surfaces  in  three-space  using  spatial  averages  of  RMS 
variations  for  surface  patches  of  width  l  measured  relative  to  a  plane-fit  to  the  local  data 
(Appendix  A).  The  average  variance  is  measured  for  geometrically  increasing  values  of  l  (/  < 
system  size  L),  and  the  dependence  of  the  variance  on  l  gives  rise  to  the  scaling  properties  of  the 
fluctuations.  The  scaling  properties  are  characterized  by  the  roughness  exponent  a  at  time  t 

wdl,t)~la  (4) 

or  is  obtained  from  the  scaling  dependence  of  the  local  height  variance,  WL(l,t),  on  l,  given 
by 


w2(/,0  =  ([/»(r,0-*/(r,t)]2)r  (5) 

where  h(r,t )  is  the  height  of  a  single-valued  surface  at  location  r  at  time  t  relative  to  a  plane-fit  to 
the  local  data,  hi(r,t)  is  the  average  local  height  at  time  t,  t  is  the  exposure  time,  and  /  is  the  size 
of  local  surface  patches.  A  DFA  is  particularly  useful  for  characterizing  the  unique  structure  of 
the  coating,  since  the  average  local  variance  of  the  substrate  can  usually  be  subtracted  from  that 
of  the  coating.  Figure  1  shows  AFM  data  and  the  patchwork  of  local  reference  surfaces 
employed  by  DFA  for  l  =  0.25, 0.5,  and  1.0-|im.  The  system  size  (scan  size)  is  L  =  5-|im. 
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/  =  0.50  /==  1.00 


Figure  1 .  DFA  reference  patches  for  a  5-|im  AFM  image  of 
sputtered  niobium.  (All  dimensions  in  |im.) 

Identifying  the  appropriate  scaling  region  in  the  Ch(r,t. )  or  wi{l,t)  data  is  critical  in 
obtaining  accurate  estimates  of  the  scaling  parameters.  The  measured  parameters  depend  on  the 
range  of  data  used  to  fit  the  linear  region.  Therefore,  we  developed  a  systematic  procedure  that 
determines  them  by  fitting  linear  and  polynomial  splines  to  the  height  correlation  results  using 
the  data  that  minimize  the  residuals  in  the  fit  (ref  9).  We  developed  this  technique  since  blind 
regression  fits  often  result  in  incorrect  values  (ref  3).  The  linear  range  is  determined  by  fitting  a 
straight  line,/;,  to  the  first  k  points  and  fitting  a  polynomial,  f2,  to  the  remaining  points  in  the  set. 
k  is  the  index  of  the  data  point  at  which  the  linear  and  curved  segments  meet.  It  is  the  variable 
knot  (index)  of  a  linear/curved  spline  and  is  a  nonlinear  parameter  determined  by  exhaustion. 

The  procedure  is  described  in  Appendix  B. 

RESULTS  AND  DISCUSSION 

A  series  of  tantalum  and  niobium  sputtered  silicon  wafers  was  prepared  for  analysis  at 
different  stages  using  a  rotating  platform  within  a  planar  magnetron-sputtering  chamber.  The 
surface  structure  of  5  niobium  coatings  (15,  30, 45,  60,  and  75  minutes  of  sputtering)  and  4 
tantalum  coatings  (15,  30,  60, and  120  minutes  of  sputtering)  was  characterized  in  terms  of 
scaling  exponents.  Each  sample  had  a  unique,  nonuniform  morphology  distribution  across  the 
surface.  Therefore,  the  region  with  the  greatest  RMS  variations  in  each  specimen  was  used  in 
the  analysis.  Simulations  of  the  physical  deposition  process  in  our  system  do  not  predict  the 
nonuniform  coating  distribution  observed  in  the  coatings  (ref  10).  However,  they  do  predict 
variation  in  the  angular  distribution  of  the  particles  arriving  at  the  surface.  Therefore,  the  coating 
structure  was  analyzed  at  different  locations  across  the  sample  to  determine  if  ctris  a  function  of 
position  due  to  particle  shadowing.  Figure  2  shows  the  variation  in  a  and  normalized  %±(60) 
from  the  center  to  the  edge  of  the  sample  for  both  sputtered  tantalum  and  niobium.  The  figure 
shows  the  wide  variation  in  RMS  surface  roughness  {%±(t)  =  V2  *  o(t))  for  both  coatings.  It  also 
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illustrates  that  although  a  decreases  monotonically  for  tantalum,  it  is  insensitive  to  position  for 
niobium.  These  results  indicate  that  particle  shadowing  is  not  a  significant  effect  and  the 
monotonic  decrease  in  a  for  tantalum  can  be  associated  with  a  phase  transition  from  /7-tantalum 
to  cr-tantalum  across  the  sample  (ref  10).  This  also  indicates  that  or-tantalum  develops  an 
inherently  more  complex  structure  than  /7-tantalum. 


Figure  2.  Roughness  exponent  (a)  and  normalized  perpendicular  correlation  length  (£j) as  a 
function  of  position  on  60-minute  sputtered  samples  of  tantalum  (Ta)  and  niobium(Nb). 

The  regions  in  each  tantalum  sample  with  the  greatest  RMS  values  correspond  to  the  a- 
phase.  Figure  3  shows  AFM  images  of  ar-tantalum  after  60  minutes  and  120  minutes  of 
sputtering.  We  were  unable  to  identify  a  unique  scaling  region  in  Ch(r,t )  for  these  coatings  (ref 
10),  as  we  did  for  the  15  to  45-minute  sputtered  tantalum  coatings.  This  is  likely  a  result  of  the 
evolution  of  the  or-coating  from  a  fine-grained  microstructure  into  a  material  with  large,  complex 
grains.  We  were  therefore  unable  to  determine  the  dynamic  scaling  exponents  for  the  evolving 
tantalum  coating. 


60  minutes  120  minutes 


Figure  3.  5-(im  scans  of  ar-tantalum  after 
60  and  120  minutes  of  sputtering. 
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The  microstructure  of  the  sputtered  niobium  evolved  differently  from  that  of  tantalum. 
Figure  4  shows  the  niobium  coating  after  t  =  15,  30,  60,  and  75  minutes  of  sputtering,  and  Figure 
5  shows  the  corresponding  C/,(r,r).  The  grain  structure  of  the  niobium  coating  is  essentially 
Euclidian  ( a>  0.88),  with  the  grain  size  (2*£//)  increasing  from  188-nm  to  368-nm.  The  values 
of  grain  size  were  validated  with  micrographs  of  the  coatings.  Figure  6  gives  the  time- 
dependence  of  £//(t)  and  %x(t),  and  shows  that  the  evolving  grain  size  and  roughness  is  consistent 
with  dynamic  scaling  theory.  The  value  of  / 3  was  determined  to  be  0.27,  and  1/z  =  0.3 1  using  all 
of  the  data  for  the  fit.  These  results  agree  well  with  published  results  for  sputtered  deposited 
coatings  (ref  11). 


15  minutes  30  minutes 


60  minutes  75  minutes 

Figure  4.  5-p.m  scans  of  niobium  after 
15,  30,  60,  and  75  minutes  of  sputtering. 
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Figure  5.  Ch(r,t)  for  sputtered  niobium 
with  t  =  15,  30, 45, 60,  and  75  minutes. 


4,(nm)  £i(nm) 

10°r — — — — — -1  r — — — 


minutes  minutes 


Figure  6.  ^//(t)  and  £jft)  for  sputtered  niobium. 

The  time  evolution  of  sputtered  niobium  was  also  analyzed  using  DFA  to  identify  long- 
range  power-law  correlations  and  correct  for  inhomogeneities  in  the  surface  structure.  We  have 
successfully  tested  DFA  on  model  sets  where  the  analytic  solution  is  available  and  used  DFA. 
We  have  also  used  DFA  in  the  analysis  of  AFM  images  where  there  is  localized  curvature  (ref 
12).  Figure  7  shows  Ch(r,t)  and  DFA  results  for  a  60-minute  sputtered  niobium  coating.  The 
DFA  results  are  scaled  to  match  those  of  Ch(r,t),  since  Irl  =  Z/2,  and  at  large  r,  Ch(r,t)  —*  2cf(t). 
The  figure  illustrates  that  DFA  produces  results  that  are  inconsistent  with  those  of  Ch(r,t).  This 
occurs  because  DFA  measures  fluctuations  about  a  plane-fit  to  local  data  and  the  niobium  grains 
(Z  <  2  *£/)  are  essentially  Euclidean  (planar  for  l «  2  *£//).  Figure  7  shows  the  breakdown  of 
DFA  for  Z  <  2  *£//>  at  approximately  300-nm.  Therefore,  the  dynamic  scaling  parameters  of  the 
niobium  coating  were  also  measured  using  fluctuation  analysis.  Fluctuation  analysis  determines 
the  scaling  properties  of  surfaces  in  three-space  using  spatial  averages  of  RMS  variations  of 
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surface  patches,  but  the  RMS  variations  are  not  measured  relative  to  a  plane-fit  to  the  local  data. 
Figure  7  shows  fluctuation  analysis  results  scaled  to  match  those  of  Ci,(r,t).  The  figure  illustrates 
that  appropriately  scaled  fluctuation  analysis  data  return  scaling  parameters  consistent  with  those 
of  Ch(r,t). 


Figure  7.  Ch(r,60),  scaled  fluctuation  analysis  (FA),  and  DFA  results 
for  a  niobium  coating  after  60  minutes  of  sputtering. 


CONCLUSIONS 

1.  The  structures  of  the  coatings  analyzed  in  this  study  are  consistent  with  those  of  a  self- 
affine  surface  fractal.  However,  the  roughness  and  correlation  lengths  cannot  always  be 
described  by  single  values  across  the  range  of  scales. 

2.  Identifying  the  appropriate  scaling  region  in  Ch(r,t)  or  )  data  is  critical  when 
measuring  subtle  changes  in  the  scaling  parameters  of  real  data.  Fitting  linear  and  polynomial 
splines  to  Ch{r,t)  or  wL(l,t)  data  provides  a  systematic  approach  for  measuring  the  scaling 
parameters. 

3.  Scaling  analysis  indicates  that  the  effect  of  particle  shadowing  on  the  intrinsic  structure 
of  the  coatings  produced  in  our  sputtering  system  is  minimal. 

4.  Detrended  fluctuation  analysis  is  a  useful  technique  for  measuring  the  scaling  properties 
of  a  complex  structure  when  there  is  inherent  bias  in  the  data.  It  is  particularly  valuable  for 
characterizing  the  unique  structure  of  a  coating,  since  the  average  local  variance  of  the  substrate 
can  usually  be  subtracted  from  that  of  the  coating.  However,  DFA  breaks  down  when  measuring 
the  scaling  properties  of  Euclidean  structures  where  variations  about  a  planar  (or  higher  order)  fit 
to  the  data  approach  zero.  In  these  cases,  fluctuation  analysis  should  be  used. 
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5.  The  c^-phase  of  tantalum  is  inherently  more  complex  than  the  /?-phase,  and  the  time 
evolution  of  its  microstructure  cannot  be  described  by  a  unique  set  of  scaling  parameters. 
However,  the  deposition  of  niobium  is  consistent  with  dynamic  scaling  theory,  and  the  measured 
parameters  provide  a  new  metric  for  quantifying  the  evolving  surface  morphology. 
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APPENDIX  A 


The  following  series  of  equations  provides  the  coefficients  of  a  plane-fit  (f(x,y)  =  a  +  bx 
+  cy)  to  the  data  Zi  that  minimize  the  square  of  the  residuals  (e,  =  a  +  bxi  +  cy,-  -  nf  for  all 
points,  m,  containing  a  local  surface  patch  of  length  /: 


m  m 


—  =  0^>a*m  +  bYJxi  +  c'£jyi=YJzi 

da  111 


TT  =  0  =>  x,  + x?  +  CZ  Wi  =  X  *.-*/ 

ob  iii i 


+»!*,  y,  tcj; 

dc  ii  ii 
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APPENDIX  B 


This  procedure  was  developed  to  determine  if  a  scaling  region  exists  and  to  determine 
estimates  of  the  scaling  region.  It  measures  the  geometric  parameters  by  fitting  linear  and 
polynomial  splines  to  Ch(r,t )  or  wL(/,r)  data,  and  estimating  the  values  of  £//and  orusing  the  data 
that  minimize  the  residuals  in  the  fit.  The  linear  range  is  determined  by  the  fitting  a  straight  line, 
fi,  to  the  first  k  points  and  fitting  a  polynomial,/2,  to  the  remaining  points  in  the  set,  k  through  m. 
k  is  the  index  of  the  data  point  at  which  the  linear  and  curved  segments  meet.  It  is  the  variable 
knot  (index)  of  a  linear/curved  spline  and  is  a  nonlinear  parameter  determined  by  exhaustion. 

The  coefficients  of  the  polynomials  are  evaluated  with  continuity  enforced  at  k,  and  the  residuals 
(£,)  are  computed  over  the  entire  range.  The  value  of  k,  for  which  the  residuals  are  a  minimum, 
is  selected  as  the  index  for  the  last  point  of  the  linear  range,  k  is  computed  using  a  quadratic, 
cubic,  and  quartic  for/2.  C1  and  C2  continuity  are  enforced  for  the  quartic  fit.  C2  continuity  is 
not  enforced  for  the  cubic  fit  to  allow  some  overshoot  in  the  data.  The  largest  of  the  three  values 
of  k  is  used  to  ensure  that  the  maximum  number  of  points  in  the  linear  region  is  selected  to 
determine  a.  If  k  is  less  than  three,  it  is  assumed  no  linear  region  exists  and  the  system  does  not 
adhere  to  fractal  theory.  £//  is  then  determined  from  the  intersection  of  the  linear  fit  and  £/.  The 
points  used  to  compute  are  at  large  r  where  no  correlation  is  observed.  The  procedure  for  the 
quartic,  cubic,  and  quadratic  fits  is  described  in  detail  in  Reference  9.  The  procedure  for  the 
quartic  fit  is  repeated  below  as 


fj(  x)  =  a+b(x-k) 

f^(x)  =  c  +  d(x-k)  +  e(x-k)^+f(x-k)^  +  g(  x  —  kf 

For  C2  continuity  at  knot  k, 

fj(k)  =  f2(k);  f'1(k)  =  f2(k);fj(k)  =  f'2(k) 
and  the  residuals  (£,)  become 

£.  =  a+bp.-y.  fori<ik 

for  i  >  ik 

=  (x.-kf 


£.=a+bp.+fq.  +  gr.-y. 

i  ri  J  6  i  Ji 


where 


p.  =  x.-k;  q.  =(x.  -k)  ;  r. 

i  i  i 


(Bl) 

(B2) 


(B3) 


(B4) 
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Minimizing  the  sum  of  the  squares  of  the  residuals  S 


S=  H  e?=Ze?=  Z  £?+  I  <=  Z<+I< 
i  =  0  ZJ?  i=0  *~*k+l  L  R 


—  =  0=>  a*m  +  bZ  P;  +  /£<?,-  +  #I>  =  £  yf- 


(B5) 

(B6) 


—  =  o  =>  a  2  p.  +  ft  Z  p.  +  /Z p.<?;-  +  #Z  p.r.  =  Z  p.y. 
to  LR  1  LR  1  R  R  LR 


~  =  0  =>  «I r.  +  *Z p.r.  +  /ZP-r.  +  gZ>f  =  Er.v. 
dg  Rl  R  11  R  11  Rl  Rl 


f|  =  0  =>  aZa.  +*>Z  p.a.  +  /Z?2  +  *1**  =  Z9,y, 

#  fl  1  Z?  1  1  Z?  1  R  11  R  1  1 


(B7) 

(B8) 

(B9) 


The  systems  of  equations  (B6)  through  (B9)  are  solved  for  coefficients  a,  b,  f,  and  g  for  k 
ranging  from  2  to  m-1.  The  data  point  corresponding  to  the  value  of  k,  for  which  residuals  are  a 
minimum,  is  used  as  the  last  point  in  the  linear  fit.  £//  is  then  determined  from  the  intersection  of 
the  linear  fit  and  %±. 
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